The giant plasticity of a quantum crystal 
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When submitted to large stresses at high temperature, usual crystals may irreversibly deform. 
This phenomenon is known as plasticity and it is due to the motion of crystal defects such as 
dislocations. We have discovered that, in the absence of impurities and in the zero temperature limit, 
helium 4 crystals present a giant plasticity that is anisotropic and reversible. Direct measurements 
on oriented single crystals show that their resistance to shear nearly vanishes in one particular 
direction because dislocations glide freely parallel to the basal planes of the hexagonal structure. 
This plasticity disappears as soon as traces of helium 3 impurities bind to the dislocations or if their 
motion is damped by collisions with thermal phonons. 
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A crystal resists elastically to a shear stress, contrary 
to a fluid, which flows. When applied to a crystal, a small 
shear stress a produces a strain e that is proportional to 
a. If the stress is released, the crystal returns elastically 
to its original shape. The elastic shear modulus of the 
crystal is the ratio \i = cr/e. However, if the applied stress 
is large, a usual crystal is "plastic" : it deforms more than 
in the elastic regime thanks to the motion of dislocations. 
With classical crystals, plastic deformation is irreversible 
and small. In metals, it requires stresses of order 10 _1 to 
10 -6 of the shear modulus and the strain rate increases 
with temperature [1, 2]. We have discovered that, in the 
case of ultrapure helium 4 single crystals around 0.1 K, 
the resistance to shear nearly vanishes in one particular 
direction while in another direction no measurable de- 
viation from normal elastic behavior occurs. This giant 
anisotropic plasticity is independent of stress down to 
extremly low values (10 -11 times the elastic shear modu- 
lus). We demonstrate that it is a consequence of disloca- 
tions gliding freely along the basal planes of the hexag- 
onal crystal structure. Some gliding along these planes 
had been observed [3, 4] in other hexagonal crystals but 
never with comparable amplitude. In 4 He crystals, we 
show that dislocations are able to move at high speed 
with negligible dissipation so that the plasticity is also 
independent of the rate at which the strain is applied. 
Moreover, this motion is reversible so that the plastic- 
ity shows up as a 50 to 80% reduction of one particular 
elastic coefficient, which we have identified as C44. The 
plasticity disappears if traces of impurities bind to dis- 
locations or if thermal phonons damp their motion. Our 
observations have been made possible by direct elastic- 
ity measurements on oriented crystals where the impu- 
rity concentration could be lowered from 0.3 ppm to zero 
at very low temperature (15 mK). The gliding of dislo- 
cations is an intensively studied phenomenon of funda- 
mental importance in Materials Science [5]. We provide 
strong experimental evidence for its precise origin. 
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FIG. 1: The shear modulus of crystal X15a as a function of 
temperature. Around 0.2 K, the shear modulus is reduced to 
72 bar, much less than the normal value (127 bar). 



We have improved some of our previous techniques [6- 
9] . Here we briefly describe them (see details in the Sup- 
plemental Material). In a transparent cell inside a dilu- 
tion refrigerator with optical access, crystals can be seen 
during growth. After determining their orientation from 
their shape (Fig. 1), we grow them in a 1 mm gap be- 
tween two piezo-electric transducers. Applying a small 
ac-voltage to one transducer (0.001 to 1 Volt, 200 to 
20,000 Hz) produces a vertical displacement (0.001 to 
1 Angstrom) and consequently a strain e of amplitude 
10" 10 to 10" 7 and a stress a = fie (10" 9 to 10" 5 bar) 
on the other transducer (/i is the relevant shear modulus 
of the crystal). The stress generates a current that is 
our signal after amplification with a lock-in. The shear 
modulus fx and the dissipation 1/Q are obtained from 
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the magnitude and phase of the signal. Various growth 
methods lead to three different types of crystals [7]: 

1 - High quality "type 1" crystals are grown by pressur- 
izing the superfluid liquid in the cell at 20 mK up to the 
crystallization pressure and staying there after the cell is 
nearly full. Random nucleation produces different orien- 
tations. In this case, there remains some liquid in cor- 
ners but not in the gap. Fresh from growth, such crystals 
are totally free of impurities, especially if one uses ultra- 
pure 4 He containing 0.4 ppb of 3 He (we also use natural 
4 He with 0.3 ppm 3 He). Slow growth at 20 mK rejects 
3 He impurities in front of the liquid-solid interface as in 
the "zone melting" purification of metals but here this 
method is extremely efficient because the temperature 
is very low [8]. However, when warming such ultrapure 
crystals, 3 He impurities come back into the crystal where 
their concentration varies with temperature. 

2 - Single crystals of lower quality (type 2) are grown at 
1.4 K. All the remaining liquid is crystallized during fur- 
ther cooling along the melting line from 26.1 bar at 1.4 
K to 25.3 bar at 1 K, below which it stays constant [10]. 
This pressure change probably creates dislocations, but 
the 3 He concentration is that of the helium gas used to 
grow them. Any type 1 crystal can be melted down to a 
small seed and regrown at 1.4 K as a type 2 crystal with 
the same orientation, or vice versa. 

3 - Polycrystals (type 3) are grown at constant volume 
starting with normal liquid at 3 K and 60 bar. 

Fig. 1 presents a measurement of the shear modulus 
H of crystal XI 5a at 9 kHz during cooling. This is an 
ultrapure type 1 crystal, although grown at 600 mK. Its 
shear modulus reaches a very small value around 0.2 K. 
At lower T, the crystal stiffens as a consequence of dislo- 
cation pinning by 3 He [6, 9, 11]. Crystals with immobile 
dislocations are stiffer than with mobile dislocations. The 
pinning temperature depends on the 3 He concentration, 
on the binding energy E3 = 0.73 K [11], and on the mea- 
surement frequency. Above 0.3 K, \i also increases, prob- 
ably because dislocations scatter thermal phonons [12]. 
We thus understand that 4 He crystals show an anoma- 
lous shear modulus between impurity pinning at very low 
T and collisions with thermal phonons at higher T. We 
call it a "giant plasticity" because it is due to the large 
motion of dislocations under very small stresses. Since it 
is reversible, the plasticity shows up as a reduction in the 
shear modulus. Note that the present measurements are 
made at kHz frequencies during hours with results that 
are independent of the sign and duration of the applied 
strain. 

Two type 2 crystals are compared in Fig. 2a: XI 5c 
(purity 0.4 ppb) and X21 (0.3 ppm 3 He), at high or low 
driving strain. When a larger drive is applied, the tem- 
perature at which 3 He bind to dislocations is lower, so 
that the plasticity domain is wider. As expected [6, 9, 11], 
pinning by 3 He occurs at lower T for ultrapure crystals. 
Fig. 2b shows the corresponding dissipation. In the plas- 
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FIG. 2: Comparison of two "type 2" crystals: X15c and X21 
grown respectively from ultrapure and from natural 4 He. The 
T-domain of the giant plasticity depends on both purity and 
strain amplitude. On graph (a), two ticks on the vertical 
axis indicate the predicted value of the shear modulus for 
immobile dislocations. The lower graph (b) shows that the 
dissipation is negligible in the high plasticity region of X15c, 
as for immobile dislocations in the zero temperature limit of 
X21. 



ticity domain of X15c, we measure zero dissipation, as in 
the low T limit of X21 where 3 He binding suppresses 
dislocation motion. The peak dissipation corresponding 
to 3 He binding is lower in temperature for the ultrapure 
crystal XI 5c and it nearly disappears at a high strain 
e = 10 _T that exceeds the threshold for 3 He unbinding 
(3 10 -8 , see [6, 13]). The dissipation above 0.3 K is inde- 
pendent of amplitude and increases with T as expected 
for a thermal phonon scattering mechanism [12]. 

We have identified the precise origin of the plasticity 
with the following quantitative study. Hexagonal crys- 
tals have 5 independent elastic coefficients [13-15] en, 
C12, C13, C33, C44 and c 6 6 = (cn-ci 2 )/2. The indices 1 to 
6 respectively refer to xx, yy, zz, yz, xz and xy, where z 
is the [0001] "c" axis. Fig. 3 shows results from 7 type 2 
crystals, all of natural purity except the ultrapure X15, 
and one polycrystal (BC2). X3 is tilted with a c-axis 
45° from the vertical. Its effective shear modulus is \i = 
0.004 C44 + 0.004 c 66 + 0.25 (cn+css-2c 13 ) , practically 
independent of C44 and cqq (see Suppl. Mat.). It shows 
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FIG. 3: Measurements (a) for different crystal orientations 
show that the plasticity is highly anisotropic. Ticks on the 
vertical axis show the predicted shear modulus with immobile 
dislocations, in excellent agreement with measurements in the 
low temperature limit where dislocations are pinned by 3 He. 



no T-dependence. Since it is very unlikely that these co- 
efficients would vary with T in just such a way to keep 
(cn+C33-2ci3) constant, we conclude that C12, cis and 
C33 are independent of T, so that X3 could be used to 
calibrate the piezo-electric coefficient of our transducers 
(0.88 Angstrom/ Volt in our first cell, 0.95 in the second 
cell). For this we used values of elastic coefficients mea- 
sured from ultrasound velocity at the melting pressure 
[14, 15] en = 405 bar, C12 = 213 bar, C13 = 105 bar, C33 
= 554 bar, C44 = 124 bar, and cqq = 96 bar. These values 
correspond to the intrinsic elasticity of 4 He crystals be- 
cause dislocations could not move at 10 MHz above 1.2 K. 
Given this calibration we could calculate the shear mod- 
ulus for all other orientations in the T=0 limit (ticks on 
the vertical axis of Fig. 3). Excellent agreement is found 
with our measurements except for XI 5 whose stiffening 
was not completed at 15 mK. For the polycrystal BC2, 
Maris' averaging method [16] leads to excellent agree- 
ment again. 

All crystals except X3 show a strong T-dependence, 
which could be due to a reduction of either C44 or cqq. 
Indeed, hexagonal crystals usually have one preferential 
gliding plane that can be either the basal plane (0001) or 
the prismatic planes perpendicular to it. For hep crystals, 
it is predicted to be the basal plane, although directional 
bonding in hexagonal metals can lead to non-close packed 
structures which favour prismatic glide [3, 4]. For X21, 
the dependence on cqq is negligible and the large varia- 
tion of C44 could be easily extracted. Assuming then that 
cqq is constant, we could calculate the variation of C44 in 
all crystals. Since X2, X5, X6 and X21 were all type 2 
crystals grown in similar conditions, one expects the T- 
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FIG. 4: Assuming that dislocations glide in the basal plane, 
one finds the same variation of C44 by 62±8 % for X2, X5, 
X6 and X21, which are crystals with different orientations 
grown in the same conditions. Assuming that they glide in 
the prismatic planes leads to absurd results (see text). 



dependence of their coefficient C44 to be the same. Fig. 4 
shows that this assumption leads to a reproducible re- 
duction of C44 (62±8 % for all type 2 crystals), although 
the shear modulus of X2 depends mostly on C44 while 
that of X5 depends mostly on cqq. On the opposite, as- 
suming that C44 is constant and cqq varies leads to absurd 
results: from 60% reduction in cqq for X5 to 300% for X6 
and to 1000 % for X2, not mentioning X21. The pref- 
erential direction for gliding has to be the basal plane. 
Only C44 depends on temperature. This is why Fig. 3 
shows a highly anisotropic plasticity. 

The amplitude dependence is also consistent with 3 He 
pinning at low T. Fig. 5 shows the shear modulus as 
a function of stress for 4 crystals at 20 mK. Above a 1 
/ibar threshold due to the breakaway from 3 He impuri- 
ties [6, 13], the shear modulus decreases. By projecting 
the stress on the basal plane, we determined the "resolved 
stress" (see Suppl. Mat.), and we obtained a reproducible 
threshold for different orientations, confirming that dis- 
locations glide in the basal planes. This threshold is hys- 
teretic because the force acting on dislocations depends 
on the pinning length Li between impurities [17]. One 
now understands why Fig. 2 shows a 3 He pinning occur- 
ing at lower T with a large strain (10 _T ) than with a 
low strain (10 -9 ). In contrast, the phonon damping is 
independent of the strain amplitude. 

Fig. 5 shows additional measurements made on X4, a 
type 1 crystal that was cooled down under high strain 
(10 -7 ) at 3 kHz in order to expell all 3 He impurities into 
the liquid. Its coefficient C44 is reduced by 80% even at 
20 mK. This value is stable in time ( 3 He stay trapped 
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ticity of 4 He crystals is mainly due to their absolute pu- 
rity at very low temperature. An interesting question is 
whether dislocations move by tunnelling or by thermal 
activation over Peierls barriers. Tunnelling is suggested 
by the reversible motion at low T without dissipation or 
strain-dependence. However, quantum fluctuations prob- 
ably reduce the kink energy by enlarging their width. 
If ksT > Ek, kinks are not relevant, as is the case for 
surface steps [24, 25] whose dynamics is also dissipation 
free in the low T limit. A calculation of could indicate 
if the dislocation motion is quantum or classical thanks 
to a negligible kink energy. 

This work was supported by grants from ERC (AdG 
247258) and from NSERC Canada. 
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FIG. 5: The relative amplitude of C44 for 4 different crystals at 
20 mK as a function of the stress projected on the basal plane 
("resolved"). At a threshold stress around 1 microbar, dislo- 
cations break away from 3 He impurities[6, 13]. This threshold 
is larger when increasing the stress than when decreasing it. 
Being free of impurities, X4 shows a reduction of C44 by 80%, 
independent of stress. 

in the liquid), and independent of the strain amplitude. 
After assuming that dislocations glide in the basal plane, 
Rojas et al. [9] found 86% reduction for a type 1 crystal 
of probably better quality (see Suppl. Mat.). Here we 
demonstrate that this assumption was correct, and that 
the plasticity is very large in particular directions. We 
have not yet prepared a crystal with zero dislocation as 
might be possible [18]. Dislocations arranged in a sim- 
ple Frank network should lead to less than 5% reduction 
in C44 [19]. In our case, they may be aligned in sub- 
boundaries [22] where they cooperate to reduce C44 much 
more. Assuming a typical dislocation density [23] A =10 4 
cm -2 , we could estimate the amplitude of their motion. 
With a 10 kHz strain e = 3 10 -8 , we find that disloca- 
tions move a distance d ~ e/Ab = 1 /im at 1 cm/s (b = 
3 Angstrom is the Burgers vector amplitude). This is 30 
million Burgers vectors per second, a giant effect with no 
equivalent in classical crystals (see Suppl. Mat.) [19-21]. 

In summary, we have demonstrated that the plastic- 
ity of 4 He crystals is reversible and anisotropic, due 
to the free gliding of dislocations parallel to the basal 
planes. Classical plasticity is irreversible because dislo- 
cations jump from one pinning site to another and mul- 
tiply, leading to "work hardening". In our case, applying 
large stresses at low T produces defects, which harden the 
crystals and are probably jogs because annealing at 1 K 
restores the high crystal quality (see Suppl. Mat.). For 
simplicity, we show measurements in this article, which 
are recorded on cooling after annealing. The giant plas- 
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THE EXPERIMENTAL CELL AND THE 
TRANSDUCER CALIBRATION 

Our experimental cell is a 5 cm 3 hexagonal hole in a 15 
mm thick copper plate, which is closed on its back and 
front faces by glass windows sealed with Indium rings 
(see Figure 1) and stainless steel clamps. It is filled with 
4 He through a thin Cu-Ni capillary (the "fill line" whose 
inner diameter is 0.4 mm). It stands at least 62 bar 
inner pressure. Inside this cell, we have glued two piezo- 
electric transducers [11] in order to shear crystals that 
are grown in the gap between them. The thickness of 
this gap was d = 1.2 mm in a first cell used for crys- 
tals from XI to X6 and d = 0.7 mm in the second cell 
shown on Figure 1, where the transducers are glued on 
their whole surface area in order to avoid any possible 




FIG. 1: The second cell used for this experiment is a 5 cm 3 
hole in a 15 mm thick Copper plate that is closed with two 
glass windows sealed with stainless steel clamps and Indium 
rings. The cell plate is attached to a dilution refrigerator 
which allows measurements down to 15 mK. 4 He crystals are 
grown inside the 0.7 mm gap between two piezo-electric trans- 
ducers by injecting liquid through a 0.4 mm capillary on the 
top (see Figure 2) . The two transducers allow applying a ver- 
tical shear and measuring the resulting stress across the thin 
crystal in the gap. A first cell was used for this experiment, 
where the gap was 1.2 mm thick. The optical access allows 
determining the crystal orientation (see Figure 3). 



bending near their edges. Figure 2 shows the "loading 
curve" , that is the stress measured as crystallization pro- 
ceeds from the bottom to the top of the transducers. 
Note that this curve shows no particular singularity at 
half loading when the liquid-solid interface passes over 
the two soldering points of grounding leads that are vis- 
ible on Figure 3. One of these soldering points is on the 
front and the other on the back so that there is no sig- 
nificant variation in the gap thickness there. Figure 3a 
shows the shape of a seed during fast growth, which is 
used to determine the crystal orientation before the crys- 
tal is regrown more slowly over the entire cell including 
the gap between the two transducers (Figure 3b). An 
ac- voltage (1 mV to IV at a frequency in the range 1 
to 20 000 Hz) is applied to one transducer, which pro- 
duces a vertical displacement u, consequently a strain 
e = u/d and finally a stress a = \ie on the other trans- 
ducer (/i is the shear modulus of the He crystal). The 
displacement u is very small - of order 1 Angstrom per 
Volt - and it needs to be accurately calibrated in order 
to obtain an absolute value for \i. The stress generates 
charges, which are collected as a current whose ampli- 
tude and phase are measured with a lock-in amplifier. 
We also use current pre- amplifiers (femto-lca-20k-200m 
and femto-lca-200-10g) as was done by Day and Beamish 
who introduced this method in their original work [2]. 
We improved their method by calibrating the transduc- 
ers' response in the following way. We first measured the 
cross talk between transducers. Although the transducer 
sides facing each other are grounded, there is a small 
cross talk between them, which mainly comes from ca- 
pacitive coupling and needs to be accurately known as 
a function of frequency. It is independent of pressure 
and temperature in our working conditions so that we 
measured it with the cell full of liquid. As explained in 
the main text, we then used a particular crystal (X3), 
which was oriented with a [0001] axis tilted by an angle 
very close to 45° from vertical. For this crystal, the re- 
sponse to a vertical shear depends mainly on the three 
elastic coefficients en, C13 and C33, with a negligible con- 
tribution from C44 and cqq. The coefficients C44 and cqq 
contribute to the shear modulus in all directions except 
for a shear at 45° from the [0001] axis. We then veri- 
fied that the measured shear modulus was independent 
of temperature in our geometry for this particular crys- 
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FIG. 2: The loading curve is the shear stress measured in the 
cell as a function of time during the growth of the crystal in 
the gap, which takes about 3 minutes in this case. Near the 
mid height of the transducers, one sees two point contacts for 
the grounding of the transducer surfaces that face each other. 
One of these contacts is on the front and the other on the back 
side of the transducer, so that there is no significant variation 
in the gap thickness there. Thanks to this design, the loading 
curve shows a linear variation with height with no singularity 
at half loading. 



tal (see Fig. 3 in the main text). Finally, we used the 
known values [7, 8] of en, cis and C33 to obtain the piezo- 
electric coefficient we needed. In the first cell it was di$ 
= 0.88 Angstrom/ Volt independently of T up to 1 K and 
0.95 in the second cell, about 5 times less than at room 
temperature. The shear modulus is fi = Id/ \ujd\ b A(j) , 
where V is the voltage applied to the first transducer, A 
is the transducers' area (1.2 cm 2 in the first cell and 1 
cm 2 in the second one), / is the current generated by the 
stress on the second transducer. In order to know /, a 
careful calibration of the gain of our amplifiers needed to 
be done as a function of frequency. In the end we ob- 
tained the absolute amplitude of the real and imaginary 
part of the response, that is the shear modulus \i and the 
dissipation 1/Q = tan((j)) where (j) is the phase delay of 
the response. For each crystal orientation, we calculated 
ji as a function of the 5 elastic coefficients (see below), 
so that we could extract the variation of C44. 



SAMPLE PREPARATION 

The quality and the purity of samples are very impor- 
tant in this experiment. The crystal quality depends on 
growth conditions as previously explained by Sasaki et al. 
[3] and by Pantalei et al. [4]. The best crystals, called 
"type 1", are grown relatively slowly (up to 50 /im/s in 



this experiment, 0.3 /im/s in the experiment by Rojas 
et al. [5]) at low T, usually around 20 mK, by pres- 
surizing superfluid 4 He up to the liquid-solid equilibrium 
pressure P eq = 25.3 bar. After nucleation on a random 
site, a crystal seed grows, falls down to the bottom part 
of the cell and the growth proceeds at constant T and P 
thanks to the mass injection into the cell through the fill 
line where helium remains liquid. This fill line is ther- 
mally anchored along its path to the cell, so that growth 
does not warm up the cell even at temperatures less than 
20 mK in this experiment. At the equilibrium or during 
slow growth, the crystal occupies the lower part of the 
cell with a horizontal surface and some capillary effects 
where it touches walls, as would a non- wetting liquid in 
a little glass bottle. This is because the growth dynam- 
ics proceeds with negligible dissipation and because the 
temperature is highly homogeneous so that there are no 
temperature gradients, only a gravity field [1]. In order 
to fill the cell with solid as much as possible, one has to 
place the orifice of the fill line at the highest point in the 
cell. For this purpose, our cells are tilted. One also has 
to avoid corners or slits where liquid would be trapped 
because of capillary effects. Finally, it is also important 
to avoid the presence of dust particles on walls because 
they are efficient pinning sites for the liquid-solid inter- 
face moving up. If one stays at P eq at the end of the 
growth, there necessarily remains some liquid in corners 
or slits, in our case at the junction between the glass 
windows and the cell body. As a consequence, the 3 He 
impurities may be trapped in this liquid if growth takes 
place at low T, because their solubility in the liquid is 
much higher than in the solid [4]. According to our ex- 
perience, it is also possible to expel all 3 He impurities in 
this adjacent liquid by applying a large ac-stress on the 
crystal at low T, because the resulting force shakes the 
dislocation and detaches them from the 3 He atoms that 
are known to travel ballistically through the crystal lat- 
tice so that they reach liquid regions in a short time and 
stay trapped there. 

We use a different growth procedure for "type 2" crys- 
tals. We grow them again at constant T and P from the 
superfluid liquid, but at 1.4 K. When the cell is as full 
of solid as possible, we block the fill line by increasing 
the pressure outside and we cool down to 1 K at con- 
stant volume. Due to the decrease of P eq from 26.1 bar 
at 1.4 K to 25.3 bar at 1 K, the rest of the liquid in the 
cell crystallizes but the quality of the final crystal, which 
occupies the whole cell now, is expected to be damaged 
by stresses. The advantage of this method is that it pro- 
duces a crystal in which the 3 He concentration is known, 
equal to the initial concentration in the gas cylinder, and 
stays at this value during temperature cycles afterwards. 
Since it is easy to melt or grow crystals by manipulating 
valves outside the refrigerator, we could melt any "type 
1" crystal to a small seed and regrow it as a "type 2" 
crystal or vice versa. 
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FIG. 3: Two photographs showing the growth shape of a seed 
crystal on the left (a), and the growth inside the gap on the 
right (b). The total width of the cell is 20 mm. The growth 
of the seed is fast enough (1 mm/s) to make facets visible on 
the left photograph (a) and to allow the determination of the 
crystal orientation (here crystal X3 with its c axis tilted at 45° 
from the vertical) . The growth inside the gap is shown on the 
right photograph (b) . Here it is our first cell where the gap is 
slightly larger, 1.2 mm instead of 0.7. The growth inside the 
gap is slower (50 /xm/s). In this first cell, the two transduc- 
ers were not glued on their entire surface and the electrical 
contacts were made on the bottom edge of the transducers. 
Between the two transducers, one can see a straight horizontal 
line, which is the solid-liquid interface outside the gap, and a 
convex one below, which shows the capillary depression of the 
meniscus inside the gap. This depression is a consequence of a 
45° contact angle with the transducer walls (for more details 
on this partial wetting, see [6]). 

"Type 3" crystals are polycrystals grown at constant 
volume from the normal liquid above the superfluid tran- 
sition. This is known as the "blocked capillary" (BC) 
method because, when cooling starts, a solid plug forms 
in the fill line near the " IK pot" of the refrigerator after 
what the growth proceeds at constant volume and con- 
stant mass (but at varying P and T of course). In this 
experiment, we started cooling down the cell around 3 
K with 60 bars everywhere. In the cell, the crystalliza- 
tion started at 2.4 K and finished at 1.7 K with a final 
pressure P= 30 bar (see "path A" in the phase diagram 
of ref. [3]). We chose these values to avoid crossing the 
reappearance of liquid near the hcp-bcc transition of the 
solid (see [3]). Indeed, recrystallization from the super- 
fluid usually ends up with a few large single crystals while 
our goal was to obtain isotropic polycrystals with small 
grains and a strong disorder. The polycrystalline nature 
of the sample is probably a consequence of multiple nu- 
cleations of seeds in a cell that is far from homogeneous 
in temperature in the absence of superfluid. We found 
that the crystals grown at the lowest temperature have 
the largest softening, that is the largest reduction of C44 
in the soft state. This is probably because their dislo- 
cation density A is smaller, with a larger free length L 
between pinning sites (a larger "pinning length"). The 
precise measurement of A and L is in progress in our 
laboratory. 

In this experiment, we reached 80% reduction of C44 
with the "type 1" crystal X4 . In a previous experiment, 



Rojas et al. [5] had found an 86% reduction of C44, after 
assuming that no other elastic constant varied, for an- 
other "type 1" crystal which was probably of even better 
quality for two reasons. First, the growth rate used by 
Rojas et al. was significantly lower (0.3 /im/s) than in 
this experiment (50 /im/s). Secondly, Rojas' cell had a 
much more open geometry with no corners and a larger 
horizontal cross section area, which allowed much more 
continuous growth without sharp jumps each time the 
liquid-solid interface detaches from some pinning site. 

In their experiment, Ruutu et al. [10] obtained small 
crystals with no screw dislocations according to their 
measurements of growth rates. Their crystals were grown 
slowly at 20 mK, but the dislocation density probably de- 
pends on the growth speed just after nucleation, which is 
difficult to monitor. We have not yet succeeded in prepar- 
ing crystals without any dislocation but it is obviously an 
exciting challenge because we expect the disappearance 
of plasticity in that case. 



ORIENTATION DEPENDENCE 

A little geometry has been necessary for the data anal- 
ysis. The elastic tensor of an hep crystal involves 5 inde- 
pendent coefficients c^- with values of i,j from 1 to 6 and 
writes: 
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The meaning of these indices from 1 to 6 is respec- 
tively xx, yy, zz, yz, xz and xy with the z-axis parallel to 
[0001], the six- fold symmetry axis - also called "c" - of 
the hexagonal structure. The coefficient cqq= (cn-ci2)/2. 
The orientation of the x-axis in the plane perpendicular 
to the c axis is arbitrary since we assume that the dislo- 
cations are distributed such that the transverse isotropy 
of the hep crystal is preserved. In our experiment, the 
axis c is tilted with respect to the vertical direction z' 
of the shear. The crystal orientation is given by the an- 
gles 6 and (j) as defined on Figure 4 where the growth 
shape is compared with a hexagonal prism. Coordinate 
transformations can be applied to quantities expressed 
in abbreviated notation using the Bond matrices [14]. 
Rotations of the coordinate system about its y-axis are 
applied with: 
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and rotations about its z-axis are applied with: 
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The elastic tensor in the transducer coordinate system 
x'yV (Figure 4) is given by: 

C 1 = M z (-0)M y (-9)CMT(-e)MT(-4>) 

where X T is the transpose of matrix X. 

The shear modulus that relates the shear strain we ap- 
ply to the component of the shear stress that we measure 
is given by: 



/i = c 44 = -(en - 2ci 3 + css)sin 2 20sin 
c^(cos 2 0cos 2 (j) + cos 2 29sin 2 (l 



cqqcos 2 (j)sin 2 6 



We used the above equation to calculate the shear 
modulus a as a function of the coefficients cij for all our 
crystals. For example we obtained: 

a = 0.0001(cn - 2ci 3 + c 33 ) + 0.933c 44 + 0.067c 66 for X2 
(0 = 89.5° and = 85°) 

a = 0.25(cn - 2ci 3 + c 33 ) + 0.004c 44 + 0.004c 66 for X3 
(0 = 45° and = 85°) 

a = 0.0001(cn - 2ci 3 + c 33 ) + 0.25c 44 + 0.56c 66 for X5 
(0 = 60° and <j> = 30°) 

a = 0.008(cn - 2c i3 + c 33 ) + 0.97c 44 for X21 
(0 = 5° and (j) = 90°) 

The values of all c^ have been obtained from ultra- 
sound velocity measurements at 10 MHz by Crepeau et 
al at 1.32K [7] and by Grey wall at 1.2K [8] c n = 405 
bar, c\2 — 213 bar, ci 3 = 105 bar, c 33 = 554 bar, c 44 = 
124 bar, and cqq — 96 bar. At such high temperatures, 
the damping of dislocation motion by thermal phonons 
being proportional to the frequency and to T 3 , [12] dis- 
locations cannot move at 10 MHz. As a consequence, 
their values correspond to the true elasticity of the lat- 
tice, without any contribution from plasticity. We have 
verified that, at low temperature in the presence of 3 He 
impurities, the shear modulus of our 4 He crystals is the 
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FIG. 4: The crystal orientation. The orientation of the [0001] 
axis "c" is denned by two angles and <j>. The axis z' is 
vertical, parallel to the transducer plane surfaces. The axis 
x' is perpendicular to the windows, and the axis y' is per- 
pendicular to the transducer surfaces. The values of and 
are obtained by matching the growth shape with a hexago- 
nal prism. In this particular case (crystal X2), the c axis is 
very close to horizontal. On the photograph, one sees that 
the crystal touches the front window so that two white lines 
correspond to intersections of the crystal with this front win- 
dow. Only the free edges of the crystal are used to determine 
the orientation. 



same as measured by Crepeau et al. [7] and by Grey- 
wall [8] . The only correction to be made for polycrystals 
grown at high pressure is the pressure dependence anal- 
ysed by H.J. Maris [9]. 

The resolved stress a r is the quantity which determines 
the force acting on dislocations. We determined it as 
follows. The strain in the transducer coordinate system 
(Figure 4) is given by: 

^=(0 e0 0) T 
In the crystal frame this becomes: 

E = M y (8)M z (4>)E' 
The stress in the crystal frame is then : 
E = CE 

The magnitude of the resultant of the shear stresses 
acting in the basal plane (E| + E5) 1 / 2 is independent of 
rotations of the crystal about its c-axis. If we assume that 
the three (1120) Burgers vectors are uniformly populated 
then we can choose the convenient orientation with one 
set of Burgers vectors along the stress resultant and the 
other two at 60° to it. Then the average resolved stress 
will be: 



c^E\/cos 2 Qcos 2 (j) + cos 2 20sin 2 (j) 

7! 
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TEMPERATURE CYCLES AND ANNEALING 
OF SAMPLES 

Most of the data that are analysed in this article have 
been obtained by cooling samples slowly (12 hours) from 
1 K down to 15 mK. The first reason for proceeding this 
way is that we confirmed that some disorder induced 
by mechanical perturbations at low T is annealed when 
warming above 0.5 K, as was previously noticed by Day 
et al. [12]. Each temperature cycle took 6 hours for 
warming up to 1 K and 12 hours for cooling down to 
15 mK, that is about one day. We used the same proce- 
dure for all crystals in order to include a recording during 
the night when perturbations from the environment were 
as small as possible (remember that we measure stresses 
down to 1 nanobar; He crystals are extremely sensitive 
to vibrations). We then transferred liquid Nitrogen every 
morning when the crystal was cold. Every three days, a 
liquid helium transfer was also necessary to keep the re- 
frigerator working. All these transfers produce mechan- 
ical vibrations that shake dislocations. After any large 
mechanical perturbation at low temperature, we observe 
some hardening of crystals, which we believe is due to 
creation of jogs on the dislocations. After annealing up 
to 1 K we found reproducible results as if jogs had been 
eliminated thanks to the diffusion of thermally activated 
vacancies. 

COMPARISON WITH A CLASSICAL CRYSTAL 

In order to compare the plasticity in helium crystals 
and in classical crystals, let us consider the historical 
measurements by Tinder and Washburn [15] on Copper. 
Qualitatively, the plasticity is the same phenomenon: it 
is due to the motion of dislocations and it is sensitive to 
the concentration of impurities. But quantitatively, there 
are striking differences. Tinder and Washburn found a 
threshold stress of about 2 g/mm 2 = 2 10 4 Pa beyond 
which a plastic strain appears in addition to the usual 
elastic response. This threshold is somewhat smaller than 
in other classical crystals, probably because of the care- 
ful growth and manipulation of these very pure samples. 
Still, it is larger by five orders of magnitude than in our 
case. In Copper, the applied stress is 0.4 10 -6 times the 
shear modulus jicu = 50 GPa. For helium crystals, we 
find a linear response (no threshold) down to 1 nanobar 
= 10 -4 Pa, which is 10 -11 times the elastic shear modu- 



lus jiHe = 12 MPa. Furthermore, Tinder and Washburn 
find a plastic strain that is 40 times smaller than the 
elastic strain so that the effective shear modulus is not 
significantly changed. In Helium, the plastic strain due 
to the dislocation motion is 4 times larger than the elas- 
tic strain, leading to an effective shear modulus that is 
reduced by 80%. In other words, the plastic response is 
2 orders of magnitude larger for stresses 5 orders of mag- 
nitude smaller than in Copper. In Copper, the plastic 
response is highly non-linear and the response time of 
order minutes at room temperature (300 K). In Helium 
at 0.1 K, the plasticity is linear so that it results in an 
effective reduction of the shear modulus, which we have 
found independent of frequency up to 16 kHz. We have 
been able to study oriented single crystals and we have 
found evidence that the gliding plane of dislocations is 
the basal plane, so that the plasticity is anisotropic. Tin- 
der and Washburn studied polycrystalline samples where 
they could not measure any orientation dependence of 
the plasticity. 
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